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Abstract 



We construct a model of a chiral transition using the well known large 
N transition in two dimensional U(N) lattice gauge theory. Restricting the 
model to a single plaquette, we introduce Grassmann variables on the corners 
of the plaquette with the natural phase factors of staggered fermions and 
couple them to the U(N) link variables. The classical theory has a continuous 
chiral symmetry which is broken at strong couplings, but is restored for weak 
couplings in the N — > oo limit. 
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I. INTRODUCTION 



Chiral transitions can occur in (non-Abelian) gauge theories coupled to fermions. These 
transitions can be continuous or discontinuous depending on the specific dynamics of the 
model. If the up and down quark masses were exactly zero, we expect the real world to 
undergo such a transition at a finite temperature. It is further believed that the order of 
chiral transitions crucially depends on the number of quark flavors |]T|. If the strange quark 
was sufficiently light, the phase transition in the real world could be first order. These 
possibilities have created a lot of interest in studying chiral transitions in lattice simulations 

Our theoretical understanding of such transitions is essentially based on universality ar- 
guments , which have worked quite well for many condensed matter systems. The present 
evidence from lattice simulations supports the predictions of universality. However, recently 
there have been some concerns j|] regarding the validity of such arguments for fermionic 
systems. Further, even if universality arguments are valid, there exist dynamics specific to 
the gauge fields that can be important near the transition and cannot be determined in the 
universal model. For example, it is well known that anomalous axial symmetry breaking 
effects can play an important role in determining the order of the transition [|l] . Recent lat- 
tice simulations || have suggested a rapid evolution in the anomalous symmetry breaking 
effects close to the chiral transition. Thus it is useful to understand the dynamics of the 
gauge fields that are responsible for chiral symmetry restoration. As a step in this direction 
it can be useful to study simplified models, that have some non-trivial gauge dynamics. 

Chiral symmetry breaking in gauge theories is closely connected to the small eigenvalue 
density of the Dirac operator. Recently some progress has been made in understanding the 
spectrum of small eigenvalues of the Dirac operator in QCD pi. The models that have 
been studied are the so called random matrix models. Such models are motivated by the 
chiral symmetries of the Dirac operator. It is possible that these models incorporate some 
dynamics of gauge fields at a crude level. When such models are extended [[jj to include finite 
temperature effects, one obtains a chiral transition. An interesting feature of the random 
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matrix models is that the dynamical variable is typically an N x N matrix. The transitions 
are obtained in the thermodynamic limit which involves taking N to infinity. In this sense, 
these can be regarded as large N transitions of a matrix model. 

A question that immediately strikes is whether there exist other known large N transitions 
that can be interpreted as a chiral transition with an appropriate extension. In this article we 
show that one can construct a chiral transition based on the large N transition discovered by 
Gross and Witten M and Wadia H] independently. They showed that the two dimensional 
U(N) lattice gauge theory with the Wilson gauge action, has a large N transition essentially 
due to the non-analytic behavior of a single plaquette integral. Consequently the pure 
U(N) lattice gauge theory with the Wilson gauge action has a third order phase transition 
separating strong and weak couplings. In the present work we show that, if we consider 
only a single plaquette, it is possible to include fermionic variables on the corners of the 
plaquette and turn the transition into a chiral transition. In the next two sections we define 
the theory more precisely and solve it in the large N limit, showing the existence of the 
chiral transition. In the final section we make some observations on the results obtained. 



II. THE ONE PLAQUETTE THEORY 

The theory we are considering is the well known one plaquette model, which consists 
of four N x N unitary matrices living on the sides of a square (plaquette). To this we add 
Grassmann N-vectors, Xn,X n associated to the corners n = 1,2,3,4. We will refer to them 
as fermionic variables as opposed to the unitary matrices that will be referred to as gauge(or 
link) variables. The unitary matrix that connects the adjacent corners nand m will be 
represented as U nm . The theory is described by the action 

S[U, X, X] = \Tr(U p + JA) - Xn( D nm + m5 nm )xm, (1) 
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where U p = U 12 U 23 U 43f U 14j . Note that the sign in front of the gauge action is natural 
from lattice QCD. The operator D couples the fermionic variables on adjacent cites. We 
will refer to D as the Dirac operator in the present model. The non-zero components of the 
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operator D nm are given by D\ 2 
with the anti-hermitian property that D 
have 

/ 



r] 12 U 12 , D 23 = V23U 23 , £> 3 4 = 7743 



t/ 43t , D 41 



(D nm Y . Thus in the matrix notation we 
\ 

(2) 

/ 

The factors rji 2 ,rj 23 ,rj^ an d Vu could be thought of as extra couplings in the theory. In 
order to fix these couplings we use the analogy with staggered fermions. Such factors occur 
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naturally in the staggered fermion formulation \TU\ as phase factors which are remnants of 
gamma matrices. In two dimension one encounters phase factors denoted by r]^(x), where 
[i = 1,2 represents the direction and x = (x±, x 2 ) the lattice site. For example one can chose 
rji(x) = 1 and r] 2 = (— 1) X1 . The well known staggered fermion action is then given by 



(3) 



In the present model since we have a plaquette, we can think of it as an object in two 
dimension. Restricting the staggered fermion action eq.(|3[) to a single plaquette the factors 
V12, V23, V43 an d ^14 can be naturally associated with rj^x). Using this connection we obtain 
one choice of the factors to be 



7712 = ??43 = VU = ~V23 = 1- 



(4) 



We will see that with the definition of the plaquette action given in eq.(|l]), the above choice 
of the 1] factorsf] is necessary for the existence of the chiral transition. 

We can now define the partition function of the theory by the usual integral over the 
gauge and fermionic variables 



[dU][d X ]W\ eM-S[U, X ,x})- 



(5) 



Other equivalent choices would lead to the same results. 
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It is obvious by construction that the theory is gauge invariant in the the usual sense. 
Further the model also has a chiral symmetry, given by 

Xn _ exp(-ie(-l) n ) X n, %„-* n exp(-i0(-l) ft ) (6) 

in the chiral limit m — > 0. This is the remnant of the chiral symmetry usually present in a 
staggered fermion formulation. This suggests that the chiral order parameter, (xx)i defined 
by 

(XX) = z J[dU][dx][dx] YsJjXnXn exp(-S[U,x,x)), (7) 

would vanish in the chiral limit. This will certainly be true for finite N. However, as we 
will show below in the large N limit, the theory breaks chiral symmetry at strong couplings. 
Further, the well known large N transition of the pure gauge theory restores the chiral 
symmetry in the weak coupling. 

III. LARGE N SOLUTION 

We will follow the work Gross and Witten M closely^ and calculate the chiral condensate 
defined in eq.(^) in the large N limit. It is important to recognize that the large N limit 
is not an approximation here, but is necessary for the transition. It is obvious that the 
above partition function will not show non-analytic properties when N is finite. Further as 
usual the chiral limit must be considered only after the large N limit is taken. Hence we are 
naturally lead to study the large N solution to the theory. 

Since the theory is gauge invariant we can go to a convenient gauge before doing the 
integrals. We will fix the gauge such that U 23 = U 1A = U A3 = 1. This naturally imposes the 
constraint that U 12 = U p . Using the properties of the Haar measure we obtain, 

(xx) = \\ [dU p ] j^Tr (7^^) exp (jpTr(U p + U}) + Tr log [D + m]) , (8a) 
Z = J [dU p ] exp (^ Tr ( U P + Up) + Tr lo § \ D + ™]) • ( 8b ) 

2 with minor differences in notation 
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The matrix D is now given by 



D 



\ 



U p -1 

-L/J 0-10 

10 1 

1 0-10 



(9) 



where it is important to remember that each element in the above 4x4 matrix is an N x iV 
matrix. The integrand in eq. fl8af ) and eq.flSbl) only depends on the eigenvalues, aij, i = 1, N 
of Up. This is obviously true for the part of the integrand that depends only on the trace 
of Up. To see that this is also true for the remaining terms that depend on the matrix D, 
note that only eigenvalues of D enter the integrand. It is easy to compute the eigenvalues 
of D. Firstly note that the only N x N matrix in D defined in eq.^), that is not a multiple 
of unity is U p . Thus using the matrix Tb n ^ m , where T is the N x N matrix such that 
[T^U p T)ij = 5ije tai ; oti G (-7T, 7r], it is possible to block diagonalize D into 4x4 matrices, 
D i , each of which is associated with one eigenvalue Q!j and is given by 



D l 



exp(zaj) — 1 

— exp(— iati) —10 

10 1 

1 0-10 



(10) 



Thus the eigenvalues of D turn out to be 



, 2(l-sin^) : A7 = ±z^/2(l-sin^), j = 1,2, 



N. 



(11) 



Note that, due to the anti-hermitian nature of D and the chiral symmetry discussed above, 
the eigenvalues are imaginary and come in pairs of opposite sign. Clearly if oti = tt = 
and D l has zero eigenvalues. This connection between a = tt and zero eigenvalues of D 
plays an important role in the existence of the chiral transition being discussed. The Haar 
measure simplifies to [dU p ] = const dT^dai A 2 (aj), where 
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J] sin 2 



(12) 



Hence the integral over T and various constants drop from expectation values such as the 
chiral condensate. Further in the large N limit, the integral is dominated by a single U p 
with a specific eigenvalue distribution. This distribution can be obtained using the saddle 
point condition, 

2 . 4sina 3 - ^ a A - a k 

— smaH — - ^ — _ , a . . = > cot . (13) 

g 2 3 m 4 + 4m 2 + 4cos 2 (^) ££2 ( J 

Except for the second term on the left hand side, this is exactly the same equation as in ||. 
The extra term comes from the extra fermionic determinant in the present model. To get 
nontrivial results we take the N — > oo limit while keeping G = Ng 2 fixed. In this limit the 
extra fermionic piece plays no role since it is suppressed by a power of N, a familiar effect 
related to the usual suppression of fermions in the large N c limit. Thus in the large N limit 
the saddle point equation, eq. (|T3D , reduces to a continuum version that was studied in |3J, 
namely 



2 r 1 

— sina(x) — P \ dy cot 

G Jo 



a(x) — a(y) 



(14) 



where a(x) is a nondecreasing function in the range < x < 1 with «j = a(i/N). This 
equation can be solved by introducing the eigenvalue density pc/ p («) — dx/da which is 
nonzero in the interval \a\ < a c , a c < ir. We have introduced the subscript U p to distinguish 
this eigenvalue density from the one to be introduced shortly, namely the eigenvalue density 
of the operator D. It is the latter that plays an interesting role in the calculation of the 
chiral condensate. Borrowing the results from || we have 
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Pu p {a) = — [1 + — cosa], \a\ < a c , G>2, (15a) 



2tt l G 
2 
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Pu ( a ) = — cos(a/2)[snT(o; c /2) - siir(a/2)]*, |a| < a c , G < 2, (15b) 
Gn 

where a c is equal to 7r when G > 2 and is given by G/2 = sin 2 (a c /2) for G < 2. 

We can use these results now to calculate the chiral condensate. In the large iV limit one 
just has to evaluate j^Tr[l/(D + m)} in the background of U p whose eigenvalue distribution 



is given by eq.(O). We obtain 
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._ , r ac . , \ { 2m 2m \ 

(XX) = da p Up {a) 2 + 2 - , (16) 

J-a. c \\ + (a) + m 2 A (a) + m 2 J 

where the values of D, namely A + (a) and X~(a) are given by eq.flTTD by just dropping the 

subscript j. 

The reason for the chiral transition is now clear. Since we need a non-zero density of 
small eigenvalues of D, we need to understand the values of a for which one obtains small 
eigenvalues of D. Note that X + (a) is always non-zero, while X~(a) can become zero when 
a = tt. Thus for the chiral condensate to be non-zero in the chiral limit, we need values 
of a close to a = it to contribute. For strong coupling, i.e. G > 2 this is clearly possible. 
However for weak coupling, i.e. G < 2, one sees that pu p (&) — for a = n, which then 
suggests the possibility of the vanishing of chiral condensate. 

We define the eigenvalue density for the operator D to be p(A). It is easy to show that 
p(X) = p(a)(\da/dX + \ + \da/dX~\), where the first term is evaluated at |A + | = A and the 
second term is evaluated at |A~| = A 

(l + |[l-2(l-A 2 /2) 2 ])/ v / 4^A2 0<A<2, G > 2, 
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' A = - 

7T 



2A 
G 



(G/2 - (1 - A 2 /2) 2 ) A!<A<A 2 , G < 2, (17) 

otherwise, 



where X x = ^2(1 - y/G/2) and A 2 = J2(l + y/G/2). In figure 1 we plot p(X) for a few 
values of G to illustrate the situation. Using the definition of p(A) we can write the chiral 
condensate to be 

<*>= f^G^O- (18) 

In the chiral limit we obtain (xx) — 77 p(0). We see that for G > 2, (xx) ~ (G ~ 2), and 
for G < 2, (%x) = 0. Further, when G = 2, we obtain (xx) ~ m - ^ we naively read off the 
critical exponents we find (3 = 1 and 5 — 1, quite different from mean field. However, since 
we do not know of any universality arguments applicable to the model, we need not take 
the exponents too seriously. Further the underlying transition is a third order transition II. 



IV. CONCLUSIONS 



Some time ago it was shown by Eguchi and Kawai (TT[ that in the large N limit, U(N) 
lattice gauge theory in d dimensions can be reduced to a model on a d dimensional hypercube. 
If this could be done in the presence of staggered fermions one would obtain a model on the 
d dimensional hypercube. In the present work we have formulated a simple model on the 
plaquette by hand. Thus it would be interesting to know if similar results can be obtained 
from the full 2-d theory in the spirit of the Eguchi-Kawai model where the final theory 
lives on the plaquette. More over we know, the underlying phase transition, responsible for 
the chiral transition just studied, is an artifact of the Wilson action used for the plaquette 
variable. The phase transition can be eliminated if one uses the heat kernel action ||12|| . This 
limits the applicability of the results obtained here. However, some qualitative features of a 
chiral transition can be motivated using this example. 

On a crude level, the results obtained here seem similar to the results of 0. First, the 
transition studied here is driven by the gauge dynamics. The fermion determinant does 
not play any role; remember that in the large N limit the fermion determinant dropped 
out of the saddle point equation eq. fllTf) . In other words, the chiral transition occurs even 
in the quenched theory. Next, the density of eigenvalues of the Dirac operator develops a 
gap for small eigenvalues above the transition (see figure 1). One might wonder if these are 
generic features of chiral transitions involving gauge fields. In QCD, it is well known that 
the number of flavors plays an important role in determining the order of the transition. 
For example, it is believed that for one flavor there is no chiral transition [|l|], due to the 
presence of the anomaly. This suggests that there will not be any transition in the quenched 
case, since the absence of the fermion determinant is likely to enhance the small eigenvalues 
of the Dirac operator and enhance the effects of chiral symmetry breaking. If this is true, 
we can expect the chiral transition to be driven mainly by the fermion determinant. It is 
through the fermion determinant that the number of flavors enters the problem. Thus the 
present model and the random matrix models do not reflect these expected features of the 
QCD chiral phase transition. 
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Another interesting observation is the role of the staggered fermion phases defined in 
eq.(f|). Notice that in the weak coupling phase the gap in p(A) for small A is related to 
the gap in the eigenvalue density pu p ( a ) at a = 7r. This can be traced back to the zero 
eigenvalues of D\ defined in eq . (|10D , which exist when q.{ = ir as discussed in section [II. 



If we had chosen all the signs in eq. (U) to be positive, the zero eigenvalues of D l would be 
related to a, = 0. Since this would not effect the saddle point equation, the density pt/ (a) 
would not change. However, since the eigenvalues of D would be related to a differently, we 
would have p(A = 0) 7^ for any G. The chiral transition would disappear. A similar effect 
can be seen by changing the sign of the gauge action. Of course we have used the natural 
sign of the Wilson gauge action that is necessary for the continuum limit to emerge from a 
lattice action. Thus the staggered fermion phases and the natural sign of the Wilson action 
in eq.(^J) have conspired together to produce the chiral transition. 

Finally, it is interesting to know if one could extend such a calculation to the more 



interesting large N phase transitions of Douglas and Kazakov |13| , which unlike the transition 
studied here is a continuum transition. Since the topology of space-time plays an interesting 
role in this continuum transition, including fermionic variables becomes non-trivial. 

V. ACKNOWLEDGMENT 

I would like to thank Andrei Matytsin for helpful conversations. I also would like to 
thank Norman Christ, Suzhou Huang and Uwe Wiese for suggestions and comments. 



9 



REFERENCES 

[1] R. Pisarski and F. Wilczek, Phys. Rev. D29 (1984) 338; 

[2] For a review of lattice results, see C. DeTar in Quark Gluon Plasma 2, edited by R. 
Hwa, World Scientific, 1995. 

[3] K. Rajagopal and F. Wilczek, Nucl. Phys. B399 (1993) 395. 

[4] A. Kocic and J. Kogut, Phys. Rev. Lett. 74 (1995) 3109 

[5] S. Chandrasekharan and N. Christ, Nucl. Phys. B(Suppl.) 47 (1996) 527-534. 

J.J.M. Verbaarschot, Phys. Rev. Lett. 72 (1994) 2531; Nucl. Phys. B426 (1994) 559. 

[7] A.D. Jackson and J.J.M. Verbaarschot, |hep-ph/ 9509324 . 



D.J. Gross and E. Witten, Phys. Rev. D21 (1980) 446. 
[9] S.R. Wadia, Phys. Lett. 93B (1980) 403. 
[10] H.S. Sharatchandra, H.J. Thun and P. Weisz, Nucl. Phys. B192 (1981) 205. 
[11] T. Eguchi and H. Kawai, Phys. Rev. Lett. 48 (1982) 1063. 
[12] P. Menotti and E. Onofri, Nucl. Phys. B190 (1981) 288. 
[13] M. Douglas and V. Kazakov, Phys. Lett. B319 (1993) 219. 



10 



FIGURES 




X 

FIG. 1. This plot of the eigenvalue density of D defined in eq.(17) for various values of the 
coupling G. At finite N, D has 4N eigenvalues. In the large N limit these can be naturally be 
described by the density given above which is normalized to 4. G = 2 is the critical coupling. 
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